AN ALGORITHMIC CHARACTERIZATION OF 
POLYNOMIAL FUNCTIONS OVER Z^n 

ASHWIN GUHA AND AMBEDKAR DUKKIPATI 
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Abstract. In this paper we consider polynomial representability of 
functions defined over Zp»i , where p is a prime and 71 is a positive integer. 
Our aim is to provide an algorithmic characterization which (i) answers 
the decision problem: to determine whether a given function over Zpn 
is polynomially representable or not, and (ii) find the polynomial if it 
is polynomial ly representable. The previous characterizations given by 
ICarlitj (ll96J) and iKempnen (|l92in are existential in nature and only 
lead to an exhaustive search method, i.e., algorithm with complexity 
exponential in size of the input. We present a new characterization 
which leads to an algorithm whose running time is linear in size of input. 
We also extend our result to the multivariate case. 
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1. Introduction 

The problem of polynomial representability of functions is central to many 
branches of mathematics. In literature, there have been attempts to repre- 
sent various functions using polynomials and power series. With the advent 
of calculus various methods were developed to approximate analytic func- 
tions using polynomials. An important milestone in this regard is the Taylor 
series, put forth by Brook Taylor. 

It is well known that if the underlying set is a finite field, every function 
from the field to itself can be represented as a polynomial. The fact that 
every function over finite fields of the for m Zp, where p i s prime, c an be 
represented by a polynomial was noted bv iHermitd ( 18631 ) . iDicksonI ( 18961 ) 
proved the above property for a general finite field. Dickson also showed 
that for a finite field of order q every function is uniquely determined by 
a polynomial of degree les s than q. Polynomials over finite fields are also 



discussed in (ICarlita. 119771) . A co mprehensive survey regarding finite fields 



can be found in ( Lidl et all 119971 ) 



In this paper we consider polynomial representability in Z^n, where p is 
a prime and n is a positive integer. Such residue rings have an elegant 
structure and their study is the first step to understand polynomial repre- 
sentability in rings. This probler n has been stu died in literatu r e and the two 
important results were given bv ICarlita ( 1964 ) and iKempnerl ( 19211 ). 

Carlita ( 1964 ) provides a necessary and suffici ent condition for a func- 



tion over Zpn to be polynomial using Taylor series. iKempneii (|l92ll ) showed 



that the only residue class rings where all functions can be represented by 
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polynomials are Zp, where p is prime. Kempner also provides a method to 
enumerate all polynomial functions over Z4 for any positive integer t. 

Keller Sz Olsonl ( 19681 ) g ive a simpler formula to express the number of 



polynomial functions in Zpn. lMullen fc StevensI ( 19841 ) provide an alternative 



formu la for the same, wh ich is also extended to polynomials in several vari- 



ables. iBrawlev fc Mullen! (|1992 ) generaliz e the formula over a Galois ring. 



Some other related work can be found in ( Zhangj . 12004 ) 



Until now the problem of polynomial represent ability has been viewed 
from a traditional standpoint and its computational aspects have been ig- 
nored. In this paper we give an alternate characterization by considering 
the set of functions over Zpn as a Zpn -module. We provide a linear time al- 
gorithm that solves the problem of polynomial representability and identify 
the polynomial which corresponds to the given function. Further, we give 
the characterization in the multivariate case. 

This paper is organized as follows. In § [2] we provide the background 
and motivation for a new characterization. §[3] contains our characterization 
of polynomial functions in Z^n . In § [5] we give an algorithm based on our 
characterization. We also discuss its correctness and complexity. In § [5] we 
determine the polynomial that corresponds to the given function. In § [6] 
we extend the characterization to functions in several variables. Concluding 
remarks are provided in § [71 



2. Background and Preliminaries 

In this section we look at polynomials over finite rings, in particular, 
polynomials over residue class rings. Let i be a positive integer. One can 
easily verify that Zt[x] is a Zj-module. Every polynomial over Z^ defines 
a function from Z^ to Z^ by the universal property of polynomial rings. 
In other words, if we allow the indeterminate x to vary over Zj, then each 
polynomial corresponds to a mapping from Z^ to Zj. Let St denote the set of 
all functions from Zj to Zf. dt — (^t)*i therefore we represent each element 
of St as a t-tuple (oq, ai,a2, ■ ■ ■ , at-i), which corresponds to the function / 
with f{i) = ai. St is a Z^-module of cardinality t*. 

It should be noted that there are infinitely many polynomials in Zf[x]. 
Let *P( denote the set of distinct functions produced by Z([x]. '^t is finite 
and a subset of St- 

Definition 2.1. A ring A is said to be polynomially complete if every 
function from A to itself can be represented as a polynomial. 

Examples of such rings are Z2, Z3, Z5. In general, Zp is polynomially com- 
plete, if p is a prime. In other words, for prime p we have *Pp = Sp- Given any 
function it is possible to construct a polynomial which corresponds to that 
function. This is achieved using Lagrange interpolation which is possible 
because Zp is a field. This does not hold for an arbitrar y integer t. Polyno- 
miall y complete structures are discussed extensively in ( Lausch fc Nobaued . 
I973I ). 



Kempneii ( 192ll ) discusses polynomials over Zj for any positive integer t. 



Kempner gives a method to compute the cardinality of *Pf The conditions 
for two polynomials to be equal as functions i.e., f{x) = g{x) mod t is 
described using the ideas of signature and characteristic of t. A method to 
enumerate all distinct polynomial functions is also provided. 

Carlita ( 19641 ) proved a key result regarding polynomial representation 



of functions in residue class ring modulo prime power. The result is very 
similar to Taylor series. The result states that a function / from Zpn to 
itself is polynomial if and only if there exist functions <l>o, *l?i, . . . , ^n-i over 
Zpn such that for all x, s G Zpu, we have 

fix + sp) = ^o{x) + (sp)$i(x) + . . . + (sp)""^$„_i(a;) mod p". (1) 



A key feature, and in a certain respect, a drawback, is that these results 
use existential proofs. The results hinge on the existence of some func- 
tions satisfying certain properties. The previous works do not address the 
issue of finding the afore mentioned functions. Consequently these results 
do not lead to any constructive method to test whether a given function 
is polynomial representable, hence the results cannot be implemented in 
computation. 

One can apply a brute-force algorithm using the result by Carlitz, by 
considering all possible functions in ^pn as shown below. 
Input: / = (ao,ai,...,apn_i). 
for all X, s G T^p^i do 

for ah $0, ^1, • • • , ^n-i e V do 

if f{x+ps) = ^o{x) + (sp)$i(x) + . . . + (sp)""^<l>„_i(x) then 
Output : / is polynomial 
exit 
else 

Output : / is not polynomial 

The above algorithm is extremely inefficient. The cardinality of the ring 
Zpn is exponential in p. \'Sp"\ = p"^ is doubly exponential in p making it 
infeasible to compute. 

We can modify the method in iKempnen ( 192ll ) to suit our problem of 



testing whether a given function / is polynomial. We can evaluate all poly- 
nomials in *Ppn and compare it with /. If / does not match any of the 
functions in *Ppn we can infer that / is not polynomially representable. The 
algorithm is presented below. 
Input: / = (ao,ai,...,apn_i). 
for all g € *Ppn do 
for all X € Zpn do 
if f{x) = g{x) then 

Output : / is polynomial 
else 

Output : / is not polynomial 
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This approach is better than the earher one, still it is very inefficient. 
^p^i, which is much smaller than 5^p", is still extremely large. One can very 
well see that given an arbitrary function in Jp" it is less likely t o be poly- 



nomia lly representable than otherwise. A simple example from iKempner 
( 192ll ) illustrates the magnitude of the sets involved. 



Example 1. Consider p = 3, n = 11. Then, 
^" = 3^1 ^ 10^ 

|^p„|= 311-3" ~ 101,000,000 _ 
|q3p„| =33-54-1076. 

We can see that even for small values of p and n, *Pp". becomes unman- 
ageably large. 

In this paper we provide an algorithm that answers the question posed 
earlier. We present a new characterization to describe polynomial functions 
over Zpn using which we bring down the complexity of the algorithm from 
doubly exponential in p to exponential in p. 



3. Characterization of Polynomial Functions over Zpn 

The question we wish to resolve can be stated as follows: 

Given a prime p and a positive integer n, and a function 
/ : Zpn — )■ Zpn, is there an algorithm to test whether / is 
polynomially representable or not? 

In order to answer the above question we make use of the module structure 
of *Pp™ . One can easily verify that for any integer t, ^t, the set of polynomial 
functions, is a Zj-submodule of ^f 

Lemma 3.1. // (ao, ai, . . . ,ai_i) and (6o, ^i, • • • , &t-i) S ^t, then (ao + 
bo,ai + bi,...,at-i + h-i) € *Pt. 

Proof. If / and g are the polynomials such that f{x) = (oq, . . . , at-i), and 
g{x) = (bo, . . . , bt-i), for X = ... t — 1, then h{x) defined as f{x) + g{x) for 
all X G Zi is also a polynomial. D 

Lemma 3.2. // (ao,ai, . . . ,at_i) E ^t, then (soq, soi, . . . , saj_i) G ^t, 
where s ^TLt. 

Proof. If f{x) = (ao, ai, . . . , at-i), for x = ... t — 1, then sf{x), which is 
also polynomial corresponds to (soq, sai, . . . , sat-i). Hence (sao, sai, . . . , sat-i) G 
*Pt. □ 

From these two lemmas we have the following proposition. 
Proposition 3.3. *Pi is a "Lt-submodule of^f 

We intend to find a 'suitable' generating set for ^p^i thereby translating 
it to a Zpn-submodule membership problem. 



3.1. Paraphernalia. 

Definition 3.4. Let f G dp"- The j cyclic shift of f, denoted by /^-'^, is 
defined as 

/<J>(i)=/(i+jmodp") 
for i = 0, . . . ,p"' — 1. 

Definition 3.5. Let vi,V2, ■ ■ ■ ,Vm £ dp"- The Zpn-submodule generated by 
Vi for i = 1,2, ... ,m and their cyclic shifts for j = 0, . . . ,p" — 1 is denoted 

by {{vi,V2,...,Vm))- 

We shall identify a set G' C *Ppn such that ^pn = {{G')). The following 
lemma helps us describe such a set. 

Lemma 3.6. //(ao, oi, . . . , ap"_i) € *Pp", then its cyclic shift (ai, . . . , apn_i, oq) € 

Proof. Let /(x) G Zpn[x] be the polynomial that gives rise to the function 
(ao, ai, . . . , Opn_i). Then /(x + 1), which is also a polynomial, gives rise to 
(ai, a2, . . . , Opn-i, ao). Hence the cyclic shift also belongs to ^p^. D 

Clearly, shifting by j places is equivalent to replacing f{x) by f{x + j). 
Hence all cyclic shifts are polynomially representable. We now state and 
prove two lemmas which are crucial in establishing our main result. 

Lemma 3.7. The function uq : Zpn — > T^pn defined as 



uo{x) = \ ' %P\'' (2) 

^ ' ^ 1 if p\x ^ ' 



belongs to *Ppi. 



Proof. We show that uq satisfies ([T]). Let ^q{x) = uq{x) and <l>j be zero 
functions for i = 1, . . . , n — 1. Now ii p \ x, p \ {x + sp) for all s € Zpn. 
Therefore, 

uq{x + sp) = 

= uo{x)+Q 

= ^o{x) + {sp)^l{x) + {spf^2{x) + ... + (sp)"-^$„_l(x). 

li p\x then uq {x) = 1 and p j (x + sp) for all s E Ijpn . 
uq{x + sp) = 1 

= uo{x)+{) 

= ^o{x) + {sp)^l{x) + (sp)2$2(x) + . . . + (sp)""^$„_l(x). 

Hence uq G *Pp'i • D 

Lemma 3.8. The function Uk ■ Zp" — > Zp" defined as 

^'=(^) = |x^ z/p!x. (^) 

belongs to ^pn for k = 1,2, . . . ,n — 1. 
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Proof. We make use of ([T|) again. Define ^q = Uk, for a fixed k £ {1,2, ... , n— 
1}. For i = 1,2, . . . ,k define $j as 



'^iix) = <, fk\k-i 



il p\ X 



For k < i < n — 1 define <^j as zero function. 

If p\ X then nfc(x + ps) = Uk{x) = and it satisfies ([1]). 

If p I X then 

Uk{x + sp) = {x + sp)'^ 

= x>^+ (^^x'-Hsp) + Qx'-\spf + ...+ (f]x'{sp)' 
= Uk{x) + (sp)$i(x) + . . . + (sp)^$fc(x) + 

= ^o(x) + (sp)$l(x) + (spf^2{x) + ... + (sp)"-^$„_l(x). 

Hence it satisfies ([1]). Therefore Uk E *Pp"- □ 



Lemma [3 .71 is in fact a special case of Lemma [3.81 when k = 0. Lemma [X 
essentially means that the following vectors can be represented as polyno- 
mials. 

uo = (1,0^_^,1,0^_^,1,...,1,0^_^) 

p-l times p-1 times p-1 times 

ui = (0,0,...,0,p,0,...,0,2p,...,(p"-p),0,...,0) 




U2 = (0,0,...,0,p^0,...,0,(2p)^...,(p"-p)^0,...,0) 

p-1 times p-1 times p-1 times 



^„_i = io,o^_^,p^~\o^_^,i2pr-\...,ip--pr~',o^_^) 

p-1 times p-1 times p-1 times 

3.2. The Characterization. We now provide the main result of this paper. 
It asserts that a function is polynomial if and only if it belongs to the 
submodule generated by Uk for A; = 0, . . . , n — 1 and their cyclic shifts. 

Theorem 3.9. / G *Ppn if and only if f G {{uq,ui, . . . ,Un-i)), where u^ for 
k = 0, . . . ,n — 1 are defined as in ^ and {{uq,ui, . . . , Un-i)) denotes the 
set generated by the vectors u^ for k = 0, . . . ,n — 1 and their cyclic shifts. 

Proof. (^^)To show that / € {{uo,ui, . . . ,Un~i)) implies / G *Pp'i. 
From Lemma 13.81 we know that Uk G *Pp" for k = 0, . . . ,n — 1. Let u^-' 
denote the j^^ cyclic shift of n^. From Lemma [3.6l we know that u^-' G ^pn 
for all A; = 0, 1, ... ,n — 1 and j = 0,1, . . . ,p — 1. From Lemmas 13.11 and 13.21 
we know that linear combinations of u^ ■^ G ^p" . 



Let / € {{uo, ■ ■ ■ ,Un-i))- Then there exist scalars a^j G Zp", for k 
0, 1, . . . , n — 1 and j = 0,1, . . . ,p — 1 such that 

/ = ao,oM^°^ + ao^iu^^^ + . . . + ao.p-i^t^^" ^ 



n— 1 p— 1 
fc=0 j=0 

Clearly all terms in the summation belong to *Ppn . Hence / € *Pp" . 

('^^)To show that / € ^pn implies / G {{uo,ui, . . . ,Un-i))- 

Let / = (ao,ai, . . . ,apn_i), where Oj G Jjp-n, for z = 0, . . . jp"- — 1. We can 

write / as 

f = VQ + Vi + ... +Vp-l, 

where Vj is the function defined as 

]ai ifi = jmodp 
I (J it z ^ J mod p, 

for J = 0, . . . ,p — 1. We now show that each Vj G ^p". From ([1]), 

aj+ps = f{j+ps) 

= Mj) + (sp)^i(j) + • • • + (sp)""^^„-i(j)- 



For j = 0, . . . ,p — 1 

f.^ /$oO') + (sp)^i(j) + --- + (sp)"^^^n-i(i) ifi^jmodp 
10 if z ^ j mod p, 

where i = j + ps. Each Vj can be written as 



V 



'J 



rii>'^ + ...+rii^l„ (5) 



where rjj: denotes the function 

(k),..^i^k{j){sp)'' iii = j + sp 
1 otherwise 

for k = 0, . . . ,n — 1. From Lemma 13.81 we can see that 

From ([5]) 

^,=cl>o(j><^> + ... + <I>„_i(i)n<i>. 

^k is well defined and ^k G Zpn for k = 0, l,...n — 1. Hence Vj is a 
linear combination of u^, for k = 0, . . . ,n — 1 and their cyclic shifts. Since 
Vj G ((uo,ui,...,M„_i)) for j = 0, ...,p- 1, / G ((uq, ui, . . . ,m„_i)). D 
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Note that not all cyclic shifts of u^ are required, for A; = 0, . . . , n — 1, but 
only the first p shifts of each Uk- This is because all the other cyclic shifts 
can be written as linear combination of the first p cyclic shifts. Hence each 
polynomial in *Ppn can be represented as a scalar sum of at most np vectors. 

This result is in fact a generalization of the generating set for vector space. 
The standard basis of the vector space 5p corresponds to uq mentioned above 
and its cyclic shifts. 



4. Algorithm based on new characterization 

Using Theorem 13.91 we provide a method in Algorithm 1 which solves 
the decision problem mentioned earlier by reducing it to a system of linear 
equations. The advantage of this reduction is that it is much easier to check if 
a system has solutions rather than check for the existence of functions which 
is done in ([1]). The linear equations can be solved by standard computational 
methods. We now present the algorithm based on the characterization. In 
the algorithm the following notations are used. 

A denotes the (n — 1) x (n — 1) matrix with elements from Z^n 



( p p^ ... p^ ^ \ 



V 


v" 


2p 


{2pr 


^p 


(3p)2 



(3p)"-i 

2 rr^ 1^^^n-l 



(6) 



\{n-l)p {{n-l)py ... {{n-l)pY^^ J 



Vi represents a p"^^-tuple which forms a subarray of input for i = 0, . . . , p—1. 
Wi represents a p"'~^-tuple of the form (p*, (2p)*, . . . , (p" — pT) for i = 
0, . . . , n — 1 i.e., Ui without the extraneous zeroes. 

One can see from ([1]) that /(x + sp) depends on f{x). In step 1 of 
Algorithm 1 we collect all the dependents in a single vector Vi of length 
pn~i_ Note that all the Vi, for i = 0, . . . ,n — l are independent of each other. 

Let q = p". Substituting s = 0, we get ^o{x) = fix) for all x G Zg. In 
step 2 we subtract this first term from each Vi to get a new vector 

(0, ai^p — Oj, aj+2p — Oj, . . . , aj+£ — Oj) 

which is written as (0, b^ ,b2 , • • • , bq_-^). 

([1]) implies that if the input function / is a polynomial then flj+ps — o-i 

(i) 

must be divisible by p. Therefore all b- s must be zeroes or multiples of p. 
With a single pass on Vi, for i = 0, . . . , n — 1, we perform this check in step 
3. If any of the Vi fails we conclude that / is not polynomial. 

In step 4, we consider the following system of linear equations over Zq, 
with variables Xj. 



Algorithm 1 Determination of Polynomial Functions 

Input: / = (ao, ai, . . . , ap™_i), where p is prime and n G N. 



Split / into p subarrays Vi such that 

^i — V'^i) Oj-j-p, (Xj-(_2p) • • • 1 0,i-\-p"—p)- 

for i = 0, 1, . . . ,p — 1 do 

tij = tij - OjlUo 

LetT;, = (0,6«,6(*\...,6(;L,)_^). 

for i = 0, 1, . . . ,p — 1 do 

for j = 0, l,...,p("-i) -1 do 
if p t ftj*^ then 

Output: / is not polynomial, 
exit 



i> Step 1 
t> Step 2 



Step 3 



for i = 0,1, . . . ,p — 1 do 



if ^ 



\ Xn-l ) 



I b? \ 



di) 



vei/ 



has no solution then 



Output: / is not polynomial. 

exit 



t> Step 4 
> A as in ([6]) 



Let $(^) = ($i*\ $^'V . . , <^]ll^) be the solution 



CO 



for i = 0,1, . . . ,p — 1 do 

n-l 

if Vi = 2~\ ^jWi then 



i> Step 5 



Output: / is polynomial. 



else 



Output: / is not polynomial. 



P 

2p 



p- 
{2pf 



p"- ^ \ ( Xl 

{2p)'^-^ X2 



I b? \ 



\ {n- l)p ((n - l)p)2 ... ((n - l)p)"-i / \ x„_i / \^ ^W ^ ^ 



di) 



(7) 
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Here we make use of ([T| to check if there exist functions ^j such that 

ai+p -ai= bf = p^i + p2$2 + . . . + p"~i$„_i 
ai+2p -ai= bf = 2p$i + (2p)2$2 + • • • + (2p)— i$„„i 



We remind ourselves that we are working with elements from the ring 
Zpn, where division by p is not defined. However, if / happens to be a 

polynomial, then all multiples of p in b- evenly cancel out. If at any stage a 
division by p is encountered it immediately implies that / is not polynomial, 
since the system has no solution. 

If solution exists for all i = 0, 1, ... ,p — 1, we then proceed to check in 
step 5 if the solution satisfies the condition for remaining components of Vi, 
i.e., we check if 

i;,G(n<'>,nf >,..., n<!.>), 

where uf^^ is the i*^ cyclic shift of Uj. If the above condition is true for 
i = 0, 1, . . . ,p — 1 we conclude that / is polynomial representable. 

The reason we choose to check for the (n — 1) components first separately 
is because had we considered all the components together we would have 
arrived at an over-defined systems of equations with p" — 1 equations for 
n — 1 variables. Computation of rank to check for solutions would take 
0((p")^) instead of O(n^) as in the case of our algorithm. 

The Algorithm 1 can be fully understood with the help of an example. 
Example 2. Consider p = 2,n = 3. Then 

Mo = (1 1 1 1 0) 
Ml = (0 2 4 6 0) 
n2 = (0 4 4 0) 

Let / over Zg be defined as 

/ = (2, 1,6, 1,2, 1,6,1). 

After Step 1: 

vo = (2 6 2 6) 
vi= (1 1 1 1) 

After Step 2: 

vo = (0 4 4) 
t;i = (0 0) 

After Step 3 we find that all the entries are divisible by 2. 

In Step 4 
For vq: 

- 2 4 W XI \ / 4 
4 M X2 j V 



for which solution exists, namely xi = 2,2:2 = 0. Hence, vq € ^Pg. We are 
misusing the notation slightly. We have avoided the extra zeroes for clarity. 
Clearly vi = (0,0,0,0) G ^s- Therefore / is a polynomial function. 

Proposition 4.1. Algorithm 1 computes whether input function is polyno- 
mially representable. 



Proof. The proof of termination of the algorithm is trivial because of the 
finite nature of the structures involved. 



From Theorem 13.91 we have that a function / is polynomial if and only 
if / G {{uo,ui, . . . ,n„_i)). More specifically, / is polynomial if and only if 
Vi G {u^^-^ , uf^-^ , . . . , u^^), for i = 0, 1, . . . ,p — 1, where u^*'* denotes the 



^th (jygi^g shift of Uj. 

In other words, there exist scalars ao,ai, 

Vi 



, Un-i in Ijpn such that 



aou'^''^ + ainf *> + . . . + an-iu'^% 



Suppose for convenience we drop the implicit zeros and write vector Vi 



lH) 



as Vi — [0^ ,02 5 • • • )''p{n-l)_;^y5 .......^ „j 

Then there exist scalars oq, ai, . . . , a„_i such that 



where b) are as described in Algorithm 1. 



/ 



b? 
bf 



Ai) 

\ V"-i)-i / 



= ao 


/ 1 \ 
1 


+ai 




yi J 


\ 



p 

2p 



\ 



( 



^n~l 



\ 



+...+a„_i 



P / 



(2p)"-i 



After step 2 of Algorithm 1, we get the first component of Vi to be zero, 
i.e., we eliminate the contribution oIuq. Let x be the vector (xi, X2, • • • , Xn-i] 
We check for solutions of 



/ 



P 
2p 



p- 



P 



n— 1 



\ 



{2py 



(2p) 



n-l 



y (p" — p) (p" — pY 



I 



Xl 
X2 



\ 



Now each Vi G ((ug 



<i> „,<«> 






/ 6« \ 



\ Xn-l ) 



di) 



9/ 



(8) 



b'-, I 
ilv / 



has a solution for x,- in 



Let A be the matrix defined in ([6]). Then, Vi G (u^*'^, n^*'*, . . . , u^^^ if and 






only if ^ • X = f j in ([7D has solution and for j = n, n + 1, 



,P 



n-l 



n-1 



Step 3 checks if A • x in ([7]) has a solution. Step 4 checks if the solution 
obtained in previous step satisfies for remaining components in ([81). D 



We now give a brief analysis of space and time complexities of the algo- 
rithm. We assume that the input is given in an array of size p", which is a 
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reasonable assumption. Also we assume that addition and scalar multipli- 
cation on vector of size p" takes 0(p") time. 

Time complexity: Step 1 takes constant time as no explicit computation is 
involved: 0(1). 

Step 2 involves a vector addition: 0(p'"~^). 
Step 3 involves one array traversal: 0(p"~^). 

Step 4 involves computing rank of (n — 1) x (n — 1) matrix to check for so- 
lution. If solution exists it can be found using Gaussian elimination: 0(71"^). 
Step 5 involves a comparison between two vectors: 0{p'^~^). 

Note that steps 2-5 can be performed in parallel as the Vi are independent 
of each other. Assuming a sequential model of computation, 

T{p,n) = 0(1) + 0(p") + 0(p") + 0(pn2) + 0(p") 
= 0(p"+pn2) 
= 0(p"+pn2). 

For all practical purposes n^ ^ p". Hence time complexity is linear in 
size of input. 

Space com,plexity: The input takes 0{p"), which is unavoidable. Apart from 
that the only space requirement is to store the (n — 1) x (n — 1) matrix which 
takes 0{'n?). Hence space complexity is 0{n'^). 

5. Determination of the polynomial 

A natural continuation of the problem is to find the polynomial which 
corresponds to the given function. This can accomplished by merely giv- 
ing the polynomials that correspond to the elements in the generating set. 
Improving upon Algorithm 1 we can obtain a solution of the system of lin- 
ear equations, if it exists. Since the solution corresponds to the scalars in 
the linear combination of the generating elements, if we are equipped with 
the polynomials corresponding to the vectors Ui defined in Lemma 13.81 for 
i = 0, . . . , n — 1, determining the polynomial of the given function becomes 
a trivial task. In this section we present the polynomials that correspond to 
the generating vectors. 

Proposition 5.1. The polynom,ial (1 — x'^'^P"') corresponds to the function 
Uf) defined in Lemma |A'. 7| as 

jo ifp\x 
[1 ifp\x, 

where (p{m) refers to Euler's totient function. 
Proof. From Euler's totient theorem we have 

for all X such that gcd(x,?Ti) = 1. 



When 771 = p" we have 

if and only if gcd{x,p'^) = 1. 

In other words we have x'^^^' = 1 mod p" if p f x for all x G Tip^. Also 
(/)(p") > n for all p > 2,n > 1. Hence {lp)'^^P"'> = mod p'', where I e Zpn, 
which means if p | x then x'^^'P ' =Q. 

From these two observations we infer that the monomial x'^'^"-' corre- 
sponds to the function 

1 if p I X. 

Then the polynomial (1 — x'^^^"'') = (p" — l)x'^'^"'* + 1 corresponds to the 
function 

11 if p I X 
which is identical to the definition of uq. □ 

It should be noted that many polynomials give rise to the function vector 
uq. The polynomial mentioned above is just one of them. It is in fact 
possible to list all the pol ynomials which correspond to uq using the method 
given bv iKempneri ( 192ll ). 



Let Uo denote the polynomial 1 — x'^^^"-*. Using uq one can easily construct 
the polynomials for all the generators of ^p". Each Ui defined in Lemma [3. 8 1 
as the function 

if p f X, 



(Xj \Jb] \ j . p I 

\ / I X it p I X 
corresponds to the polynomial Uj given as follows. 

,• fo ifpfx 

I X II p I X 

The cyclic shifts of Ui are obtained by replacing x by x + j in each Uj. 
The polynomials corresponding to the generators are 



UQ = l- x^(P") (9) 

ni = x^(l-x'^(P")) (10) 

np^(x + j)(l-(x + j)<^(P")) (11) 
for i = 1, . . . ,n — 1 and j = 1, . . . ,p — 1. 

Written explicitly the desired polynomials are 

1 - x<^(P"), 1 - (x + 1)'^(P"), . . . , 1 - (x + p - 1)'^(P"), 

x(l-x'^(P")),(x + l)(l-(x + l)<^(p")),...,(x+p-l)(l-(x+p-l)'^(p")),..., 
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^)). 



6. Polynomials in several variables 



The problem of determining whether a given function is polynomial can 
be extended to functions over several variables as well, i.e., given a function 
/ : (Zpn)"^ — y Zpn, where m is a positive integer, can we determine whether 
/ can be written as a polynomial in m variables? The ch a racter ization given 
in ([1]) is extended to multivariate functions in ICarlita ( 19641 ) . As in the 



case of single variable the ch aracterization is existent ial in nature. Some 
related work can be found in lKaiser &: Nobaued ( 19871 ). We show that our 



characterization Theorem 13.91 can be extended to multivariate functions. 



Let ^j^ denote the set of all functions from (Zpn)"^ to Z^n. Let *Pp™ 
denote those functions which are polynomially representable. 

The definition of cyclic shift in mulitvariate case is non-trivial, but follows 
closely the univariate case given in Definition 13.41 

Definition 6.1. Let f G ^pn . For (ji,..., jm) S Z*" we denote its cyclic 
shift by J<Ji'---j™> ayj(f define it as 

/^^^'■■■'^"'^{xi,...,Xm) = fixi+ji modp'^,...,Xm+jm mod p^) 
for all (xi, . . . ,x.m) G (Zpu)"^. 



For functions involving several variables the result bv lCarlita ( 19641 ) given 



in ([T]) takes the form: / : (Zpn)™ — y Zpn is polynomial if and only if there 
exists suitable functions <5ii,...,i„ : (Zpn)"^ — y Zpu such that 

f{xi+pSi, . . . ,Xm+pSm) = ^ $ii,...,i„ (xi, . . . , Xm)(pSi)'l . . . (pSm)*" mod p". 

ii+...+im.<n 

Using the above result we define a generating set similar to the one defined 
earlier in Lemma [ 



Lemma 6.2. The function Uki,...,km '■ (^p")™" — ^ ^p" defined as 

x^...x^ if p\ Xi for all i = 1, . . . ,m 



Uki....,km [Xi, . . . , Xfn) — In -^l.^ j. i ^ -i 

[ (J if p\ Xi for at Least one i = 1, . . . ,m 

belongs to ^^ , where < ki, . . . ,km < n. 



Proof. Proof is by induction on m. For m=l, the above statement is true 
by Lemma 13.81 Assume that the function Ufcj^,,.^fc^_-^ : (Zpn)™"^ — > Zpn 
defined as 

( \ — J ^1^ • • • ^m-i if p I Xj for all i = 1, . . . , m — 1. 

Uk,,...,k^-Axi, ■ ■ ■,xm-i) - I Q otherwise, 

is polynomially representable. Let h G Zpn[xi, . . . ,Xm-i] be the polynomial 
which when evaluated over (Zpn)"*"^ gives the function Wfc^^...^fc^_j. 



Consider the function u'^ : 'Lp-n 



\n defined as 



wi (x) 



ifpfx, 



c " if p I X. 

From Lemma 13.81 we know that it is polynomially representable. Let g S 
TLyn \xrr\ be the polynomial that corresponds to the function u'^ . 

Consider h and g as polynomials in Zpn [xi, . . . , Xm] . Clearly hg £ Zpn [xi, . . . , Xm] 
Let / = gh. As a function / is defined as 



Xi . . . X 



f{xi,...,Xm) = < 



™— 1 „fc„ 



if p I Xj for alH = 1, . . . , rn, — 1 and p | Xm 
x^ . . . x^Z\ -0 if p I Xj for alH = 1, . . . , rn, — 1 and p \ Xm 

0-0 



if p I Xj for some i = 1, . . . , ttz — 1 and p \ x„ 
if p I Xj for some i = 1, . . . ,m — 1 and p\ Xr, 



That is we have 

/ i^l) • • • 5 ^m j 



J x^^ . . . x^ if p I Xj for all i = 1, . . . , m, 
1 if p f Xi for some i = 1, . . . ,m, 

which is identical to nfci,...,A:„(a;i, . . . ,Xm)- Hence nfei,...,fc^(xi, . . . , x^) G 

Theorem 6.3. / G ^3^"^^ i/ and only if f £ ((■Ufci,...,^™ : fci + . . . + fem < n)), 
where {ki, . . . , km) G Z" aJ^f^ Uki,...,k„i is defined as above. 

Proof. The proof is similar to one given in Theorem 13.91 One implication is 
trivial. To prove that / is polynomial implies / G {{uk-^^...^km, ■ ^i + - • • + ^m < 
n)) write / as 

o<ii,. ■■■im<P 

such that fji,...jVn • (^p")*" — ^ ^p" defined as 



■^JlvJ. 



^(ai, . . . ,am) 



/(ai, . . . , am) if Oj = ji mod p for all i = 1, 
otherwise. 



,m 



That is /(ai, . . . , am) is placed in exactly one of the p™ different i'ji,...jVn 
functions. We now show that each Vji,...j„ G ^p™ for all (ji, . . . ,jm) G ^™'- 

This can be written as 

m 
^ ^ki,...,k^{ji,---,Jm)Yl{psi)''' ifai=ji+psi 

ki+...+km.<n i=X 

for alH = 1 , . . . , 77Z 
otherwise. 



Vh,...,j^{ai,...,am) 



Let riku-^k^ '■ (^p")^ 



■^p 



,n such that 



u 



'Jlr--j™ ~ 2^ 'nki,...,kml 

ki+...+km<n 
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where 



%i,...,fc™(ai,--- ,am) = < 

for all (ai, . . . , am) € (Zpn)™. We have 



^fci,...,fc„(ii, • • • , im)(psi)^i . . . (ps^)^™ if ai = ji +psi 

for alH = 1, . . . ,m 
otherwise. 



11ki,...,km — ^fci,...,A:mUl' • • • ) Jmj ' '^k\,.].,km ' 

where 'Ufc"'^'"^"''" denotes the (ji,...,jm,) cyclic shift of Ufci,...,^™- Hence 
each tiki,..., km € {u'^ll]]]^^^ : ki + . . . + km < n). In other words, ??fci,...,fc„ e 
{{uki,...,k„i '■ ki + . . . + km < "-)), which implies fji,...jv„ and therefore / € 

((iffci,...,fc„ : fci + ... + fem < "-))• □ 

Using the above result we can obtain an algorithm similar to Algorithm 
1 that determines whether the given function is polynomial or not. The 
complexity in multivariate case is 0((np)™), which is linear in the size of 
the input. 

Determination of the polynomial is extended to the multivariate case in a 
natural way. In the case of m variables we know that ^l^ is generated by 
{uki,...,km\ ki + . . . + km,n}. The function uo,...,o is given by the polynomial 

(1 -xf-"))(l-xf "))... (l-x^^). 
In general the function vector 

{x^^ . . . x^ if p I Xj for alH = 1, . . . , m 
otherwise, 

is given by the polynomial 

X^ ^2 . . . U^m \^ -^1 )\^ -^2 ) ■ ■ ■ \^ -^m )■ 

This way it is possible to determine the polynomial that corresponds to the 
function in multivariate case as well. 

7. Concluding Remarks 

In this paper we have considered the problem of polynomial represent abil- 
ity of functions over Z^n. A new characterization of polynomial functions 
is given that leads to a non-exhaustive algorithm which runs in linear time. 
We have also given a method to identify the polynomial that corresponds to 
the given function by providing the polynomials for the generating vectors. 
The results are extended to multivariate case as well. 
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